epidataBL[, region : = federal_countryinfo[match_hk (epidataBL$ CountryExp, federal_countryinfo$ dataname), region]]
epidataBL[, NewConfCases7_100k7 : = 7 * NewConfCases7/ Ntotal* 100000 ]
epidataBL[, AllDeaths_100k : = AllDeaths/ Ntotal* 100000 ]
epidataBL[, is_Germany : = grepl ("Germany" , region)]
epidataBL[, Altersgruppe2 : = ifelse (Altersgruppe== "all" , "all age groups" , Altersgruppe)]
p_epi = ggplot (epidataBL[Altersgruppe == "all" & DateRep<= maxdate],aes (DateRep, ifelse (NewConfCases7_100k7< 1 , 1 , NewConfCases7_100k7), color = region) ) +
ylab ("7-day-testpositives \n per 100 000" )+
geom_line () +
xlab ("" )+
labs (color = "" )+
# facet_wrap(~Altersgruppe2 ,scales = "free") +
scale_y_log10 (breaks = log_breaks (10 ),label= label_comma (accuracy = 1 ) ) +
theme_minimal (base_size = 14 ) +
ggtitle ("A)" )+
# scale_linewidth_manual(values = c(0.7,2))+
scale_color_manual (values = col_vector)+
scale_x_date (labels = date_format ("%b-%y" ))+
xlab ("" ) +
theme (legend.position = "top" ,
strip.text.y = element_text (angle = 0 ),
axis.text.x = element_text (angle = 45 , hjust = 1 , vjust= 1 )
)
p_epi